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V 

RESEARCH Ma^ORANDUM 

TEE WAVE DRAG OF ARBITRARY CONFIGURATIONS IN 
LINEARIZED FLOW AS DETERMINED BY AREAS 
AND FORCES IN OBLIQUE PLANES 
By Hsirvai^i Lomax 

SUMMARY 


The ■wave drag, based on linearized theory, of any lifting or non- 
lifting planar or nonplanar object in a steady supersonic flow is shown 
to be identical at a fixed Mach number to the average wave drag of a 
series of equivalent bodies of re'volution. The streamwise gradient of 
area, measured normal to the free stream, at a section of each of these 
equivalent bodies of revolution is given by the sum of twa quantities: 

1. The streamwise gradient of area, measured in parallel oblique 

planes "fcangent to the Mach cones, along the given object 

2. A term proportional to the resultant force on the object measured 

in the same oblique planes. 


INTRODUCTION 


The evaluation of the wave drag of airplanes traveling at transonic 
speeds has, recently, been greatly simplified by the discove37y of what 
has been termed the ’transonic area rule which ean be stated as follows: 

A (reasonably smooth) airplane flying at a 
near sonic speed has the same drag as a body 
of revolution v±th the same cross-sectional 
area in planes normal to the flight direction. 

The original statement of this rule (ref. l) has since had many eacperl- 
mental verifications. From a theoretical ■viewpoint, a paper written by 
Hayes (ref. 2) also contains the concept in the sense that it is a result 
of linearized supersonic flow theory when the Mach n-umber is allowed to 
approach 1, More recently, several papers have developed theoretical 



2 


NACA RM A55A18 


techniques which make it possible to include the study 9 f thickness 
ef,fects in slender airplane theojiy^ a linearized theory devoted 

specifically either to airplaneh flying at any speed but geometrically 
slender in the streamwise direction, or to airplanes flying near the 
speed of sound but smooth enough and having an aspect ratio low enough 
for the pressure coefficient to be small almost everywhere on their 
surface. The rule stated above lies at least implicitly in published 
restilts based on the latter theory. See, for example, references 3 4. 

When linesirlzed supersonic flow theory was shown to be consistent 
with the transonic area rule as the Manh nimiber approached unity, the 
question naturally arose as to whether or not it could be used to develop 
an analogous rule applicable to the evaluation of wave drags on airplanes 
flying at supersonic “bpeeds. Thereafter, a formula, often referred to as 
the supersonic area rule, was preseiited" (ref s . 5 and 6) which did extend 
the transonic area rule, again expressing the drag in terms of the stream- 
wise distribution of airplane cross-sectional area - but this time in 
terms of cross sections taken in planes tangent to the characteristic 
Mach cones. However, this supersonic area rule, when applied to general 
airplane shapes, gives only an approximation to the correct llnearized- 
theoiy value of the wave drag. The degree of approximation ranges from 
very good, for cases such as a nonlifting wing centrally mounted on a 
slender body (essentially the case to which Jones (ref. 5) limited his 
result), to very poor, for cases such as the .Busemann biplane or a shro\ided 
body similar to that studied by Ferri in reference 7« Further, the quan- 
titative error made in applying the rule to general shapes flying at 
supersonic speeds was unknown. 

The object of this report is to piresent a formula which gives the 
correct (linearized theory) value of the wave drag for any lifting or 
nonlifting object in a steady supersonic flow as the average of a series 
of wave drags of equivalent bodies of revolution 


LIST OF IMPORTANT SYMBOLS 


D 

l{x,9) 

M 

n 


local pressure - free -stream pressure 
pressure coefficient, * * 

q 

wave drag 

oblique section lift (see eq. ( 17 )) 
free- stream Mach nuinber 
normal to airplane surface 
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Si 
Sa ‘ 
S 3 

• 

q 

Tl 

s(x,e) 

Uo 

2=0 

x,y,z 

x,r,0 

P 

V 

i,cr 

*^0 

Po 

9 


direction coslneB with respect to the x,y,z axes 

1 2 

free- stream dynamic pressure, — PqUq 

s/Cy-yx)^ + (x-xi)® 
oblique section aj?ea (see eq. (l2) ) 
speed of free stream 
X - pr 

Cartesian coordinates in wind-axes system 
cylindrical coordinates in wind-axes system 
n/m^ - 1 

conormal to airplane surface 
coordinates defined in eqiiatlon (lO) 
coordinate defined in sketch (d) 
free- stream density 
perturbation velocity potential 


DEVELOPMENT 


Consider an airplane in a steady supersonic stream. Let the surfaces 
of the airplane be Inclined to the free stream at angles small enough for 
the dlstiirbed flow field to be adequately approximated by solutions to 
the wave equation 


P^xx “ ^Pyy ~ Tzz ^ 

where the free stream is moving in the positive x direction and 

= M^ - 1. A general solution to equation (l), developed by Volte27ra, 
is given by 


cp(x,y,z) 



xx +7(x - xx)^ - p^i^ 
^rx 


dSi 

( 2 ) 
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where ri^ = (y - Yi)^ + (z - Zj.)^, dS^ is an element of surface area 
on the airplane, Vi is the outward conormal (the conormal to the char- 
acteristic cone lies along the cone) to that element, and t is that 
portion of the airplane surface 5d.thin the Mach forecone from the point 




Sketch (a) 


We wish to use eq\iation (2) to cal- 
culate the wave drag of a given object. 

For this purpose the wave drag is esipressed 
in terms of the perturbation velocities 
induced by the object on €ua enclosing 
cylindrical control surface of infinite 
radius. This equation (see, e.g., ref. 8) 
can be written in terms of the cylindri- 
cal coordinate system defined in 
sketch (a) (notice that the control sur- 
face is parallel to the free-stream 
direction, that is, we are using wind _ 
axes) as 


D = -Po f de r dixf lim cpr 9 x') (3) 
'-'o *^-00 \r— :»oo / 


The wave drag of an arbitrarily 
shaped object (arbitrary except, of course, 
for the condition that equation (l) must 
adequately represent the flow field) can 
be calculated by inserting equation (2) 
into equation ( 3 ) . 


First let us inspect the conormal derivative 8/8vi and the differ- 
ential area dSx- By definition (see ref. P« 2h6) 


b 


2- 8 _ 8 _ 

-P'Tai ^ — + na — + ng 
axi 8yx 





+ 



where rix, Zs, and Hg are the direction cosines of an outward normal to 
the airplane surface . For surfaces making small angles with the free 
stream 


ng^ + *» 1 

Sx « streamwise surface slope 
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Hence, within the framework of linearized theory 


Svi 

where nj, is either the normal to 
the E-urface in the yi,Zi plane. 
dSi can he es^ressed as 


on.2. 

the airplane surface or the normal to 
Similarly, the differential area 


dsidxi 

dSx = . ■ r 

^1-5^ 


where dsi is an element of arc along the airplane surface in an 

Xi = constant plane; so, again, within the framework of linearized theory, 

one can write 


dSx = dsxdxi 


(5) 


For convenience equation (2) is divided into two parts such that 


cp(x,y,z) = qpx(x,y,z) + cp 2 (x,y,z) 


( 6 ) 


where 


9i(x,y,z) = - i Ay' ycpn^(xi,Bi)ln 


7(x - Xi)2 - 




dxxdsi 

(7) 


and 


<P2(x,y,z) = f / cp(xx,sx) la 

2k oxO^J onx 


X - Xi + 


7(x - Xi)2 - 




dxidsx 

( 8 ) 


Consider, first, equation (?)• According to equation (3), we need 
to find the vriue of <Px and in the limit as r (which equals 

vy® + z^) goes to infinity. However, since no disturbances can he 
induced ahead of the foremost Mach cone enveloping the disturbing object, 
it is convenient to increase x as r is Increased so that the point 
x,r,0 remains in the vicinity of this Mac^ cone. Therefore, we first 
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set (see sketch (b)) 

X = Xq + Pr 

and then let r approach Infinity. 
In this way, as r becomes very 
large, one can show that equation ( 7 ) 
reduces to 


Sketch (b) 


cpi(x,r,e) 


-1 

2n 2pr 



<Pai(xi,Si)dxidsi 

n/xq - Xi + py^cos e + pzisin 9 


(9) 


By means of the transformatione 


I = Xi - pyxcos e - pzisin 0 

O’ = Si 

equation (9) further simplifies to 


( 10 ) 


_ o Xq 3 1 r> 

cpi(x,r,e) ■-—= / , / cpni[£ - f(g ,cr) ,ff] do (ll) 

2jt v/2pr'J-oo n/xo - ^v/qc 


where / da is a line integral, around the airplane surface in the 
'^oc 
oblique cut. 


The velocity potential q)i given by equation (ll) is exactly the 
same as that Induced on a large cylinder by a line of sources distribu- 
ted along the Xi axis from -00 to Xq , the variation of their strengths 


being given by / cpmda. (This was first pointed out by Hayes in 

WQC p 

ref. 2 .) The physical significance of the term / cpn^da is 

One 


more 
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readily grasped by referring to sleet cb (c) . Imagine a series of Mach 



planes parallel to the yx axis each given by the equation 
Xi - pzx = constant. Place the airplane in its normal flight attitude. 
Each Mach plane slices through the airplane, defining, thereby, a certain 
area composed of the region on the Mach plane within the airplane surface. 
Project these areas on planes normal to the free stream (l.e., yxj^x 

planes) and designate the resulting area distribution by S(|, ^), The 

integral / <pnx<^ then proportional to the streamwise rate of change 
'^oc 

of these normally projected, obliquely cut areas; that is, for the air- 
plane so placed. 


= uo A s(e, I) 


Now, keeping the Mach planes fixed, ^ revolve the airplane about the Xx 
axis (not about Its own body axis unless the latter happens to coincide 
with the Xx axis in our wind-axes system) and repeat the above process 
for all orientations in a complete 360 ° rotation. For any given angle 


= Ho ^ S(|,0) = UoS*(|, 0 ) 


( 12 ) 


^Or holding the airplane fixed, rotate the ifach planes - always, 
of course, keeping them tangent to the characteristic ifach cone. 
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Since one can show combining equations (ll) 

and ( 12 ) yields 

(,3) 

r-3»oo P r^oo 2it 'v/2pr - | 


and, by means of equation (3)^ this gives the complete contribution to 
the wave drag of the first term in equation (2) . 

Consider next equation (8) . Taking the derivative with respect to 
ni as indicated, we have 

, ,/dyi dzi \ 

<p(xi,si)(^^ cos e + Bln ^dxidsi 

[(y-yi)^ + (z-zi)^l 3^(y-yi)^- P^(z-zi)^ 




Proceeding as before, setting x=Xo+Pr, y=r cos 6, z = r sin 0, 
and letting r go to infinity, we find 


cp2(x,r,0) = 


P 


2 n •J 2 ^T 8 xo 


cp(xi,Sx)(^ cos 0 + ^ sin ^dxidsi 

sj Xo - Xi + pyiCOB 0 + pz^sln 0 


(14) 


Introduce the transformation given by equation (lO) and q >2 can be 
expressed as 


cp2(x,r,0) = 


p 8 


. ... r <p[|.f(g^(j)^(j/^coe 0 + ^Bln ^do 


2n */^r 8x0 J -03 Vxq-I '^oc 


Notice (see sketch (d)) that the term ® ■*" ^ds^ 

is simply the component of dsi normal to the plane 0 = constant. 
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I 


Designate this direction by the 
coordinate Og , as in sketch (d), 
and equation (l4) becomes 


qpa(x,r,e) 


„ p ■ --A r ^ r 

2a -s/2|3r *7x0-1 ^oc 



Integrating by parts (notice that 
— q)[| - f(g,d),cr] = u[| - f(g,cr),ff] 

where u(x,d) = — cp(x,tf), since in 
ox 



dz, 


linearized theory csn be COS/‘GOS0 -f- Siti^SinS ^ COS(^-^) 

neglected relative to unity) and using 
the relation for pressure coefficient 



Ihe velocity potential <pa is again exactly the same as that induced 
on a large cylinder by a line of sources distributed along the Xj_ axis 
from. -00 to Xo, the variation of their strength this time being given by 



The physical significance of 




CpddQ can be demonstrated 
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with, the aid of sketch (e). As before, imagine a series of planes 



(xi - pZjL = constant) parallel to the eocls. Again place the aiiTplane 

in its normal flight attitude. TPhen if we define l(|, it/2) as the lift 
(the ccmiponent of net resulteint force in the Zi direction, positive 
upward) on a given section formed by the intersection of a Mach plane 
with the airplane svirface, one can show, for the airplane so placed. 


where q is the free-stream. dynamic pressure. If we keep the coordinate 
system and Mach planes fixed and revolve the airplane about the Xx axis. 


for each 9 the term q 


/ Cpdcg 
*^00 


lift^ on the obliquely cut section 


represents at a given 
. In general. 


X the net 



I Ul,e) 


(17) 


Now combining equations (l6) and (iT)jt one can show 

®If the airplane is fixed and the Mach planes are rotated, l(i,9) 
represents the resultant obliquely cut section force normal to the free 
stream and parallel to the plane 9 = constant, and dag is an element 
of length normal to that plane. 
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(cpa^) 


r-a,oo 



r-»(» 


Up 

at s/2pr 





( 18 ) 


where i*(|,0) = — that is, the streamwlse gradient of the lift 
^6 

on the obliquely cut section. Finaily, if equations ( 13 ) and ( 18 ) are 
placed into equation ( 3 ) and the x integration is carried out, there 
results 



where for constant 0 the intersecting Mach planes extend from 
X = -Li(0) to X = L(e) . 

Equation (l9) gives the wave drag of any lifting or nonlifting air- 
plane in a steady supersonic stream, the only approximations being those 
basic to linearized theory. 


SOME SIMPLE APPLICATIONS 
A Plane Wing 

In order to demonstrate the general applicability of equation (l9), 
let us use it to solve some sinqjle problems. Some of the simplest kinds 
of aerodynamic problems are those concerning single planar systems, that 
is, systems composed of thin wings, the sirrfaces of which are everywhere 
close to a given plane. For such cases the area and lift terms in equa- 
tion ( 19 ) have certain symmetry properties, ^eciflcally 


s*«(x,0) = s**(x,-e) 

z*(x,0) = - 2*(x,-0) 
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or i'(x,e) Is an odd and S’*(x,0) Is an even fimction of 0. Hence, 

pSJt 

J S* ' (x,0)Z*(x,0)d0 = 0 

and equation ( 19 ) reduces to 
_ , r>zsc p^(Q) pL(0) 

£= -^ / d0 / dxi / dxaS»'(xi,0)S«*(x2,0)lntxi-X2l - 


k:x- 


pL(e) 

pL(0) . 



^-Li(0) ^ 

-Li(e) ^ 


Z*(xi,e)Z* (X2,e)lnixi-X2( (20) 


which simply exemplifies the well-known: fact that in the study of a thin 
planar wing, lifting and tlilckneas affects can he analyzed separately. 
Further, since one can prove (for any f(x) that can he e:qjanded in a 
Fourier series) 

f f ( Xi) f ( X2) In I X 1 -X 2 1 dxidxa > 0 



eqiiatlon (20) shows that the only way the wave drag of a single wing can 
he zero is for hoth S**(x,&) and 2*(x,0) to be Identically zero, which 
is the trivial case of no wing at all. 


A Biplane 


Consider next a biplane, in partic\ilar, the Busemann type biplane 
shown in sketch (f) traveling at the hipest Mach number at which the 
wave drag is kncfwn (a priori) to be zero, that 1 b, at a Mach number 


equal to 



Since the flow is two-dimensional, we need only 


study the cut represented by 9 - ^ and consider strips of unit width 
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in the spartwlse direction. Sketch (f) 
shows the variation of S(x) , S*(x), and 
/ % 

■PT~ tx) with X. One can show that the 
2q 

two terms S*(x) and ^(x) are exactly 

equal for any given value of x and 
thus_, hy equation ( 19 )# the total wave 
dra^ Is zero - the correct result 
according to linearized theory. 


s'(x) 

S(x) 






~cz 


tz 



Xo 



L_ 


L 








Shrouded Bodies of Revolution 


As a final example, let us develop 
a method for calculating the shape of a 
hody of revolution which when enshrouded 
hy a thin cylindrical tube (having neg- 
ligible thickness) will result in zero 
wave drag for the combination. Let the 
length of the tube be c and its 
radius R (see sketch (g)). 

Assume a certain loading on the 
cylinder (l.e., a value of 
ACp(x) = - Cpouter^* 

equation of the projection in the —r 

xy plane of the curve formed by | 
the intersection of the oblique 
plane, x = | + pz, and the cylin- 


h . 




/^Mach 

lines 


2R 


Sketch (f) 

— Mach cones 
pylindn'cal shroud 
Body of 

revolution 



der, y^ + z® = R^, is given by 










Sketch (g) 


Hence, the oblique section lift can be written 


UP 

q 


= 2 / ACpU + p Vr 2 - y2)dy; -pR < S < 0 

'Jo 


(21a) 
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= 2 r ACpCi + p - y2)dy - 2 f ACp(| - p n/r^ - y2)dy; 

O / AA2 


n 


0 < I < pR (21 Td) 


= 2 r [ACp(| + p n/r^ - 7 ^)dy - ACp(g - p 7 r® - 3^)1 dyj 
J Q 

PR < I < c - pR (21c) 


and so forth. Now If the Interior body of revolution is slender enough 
to be analyzed by slender-body theoiy (i.e., source strength Is propor- 
tional to S’(x)), the oblique section loading on the body Itself can be 
neglected and the body shape for zero total wave drag is, therefore, 
simply 


s’(i) = ^ lii) 

2q 

where the values of l(|) are given by equations (2l) . 









Sketch (h) 


If the body determined by the above 
process is real and closes, the enshroud- 
ing tube will be a perfect cylinder and 
the velocity of the air everywhere exter- 
ior to the shaded region in sketch (h) 
will be eqTial in magnitude and direction 
to the free- stream velocity,® 

The solution given by equations (21) 
and (22) has been applied to a tube for 
which 


AOp = 


(f)(-s) 


®Thls follows directly from a theorem due to E, W, Graham (unpub- 
lished) which states that the flow fields of all finite systems having 
a minimum wave drag under the same specified restraint are identical out- 
side their enclosing fore-and-after Mach cones. 
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wiiere A is a constant. There resttlts 


S‘(x) 


S»(x) 



and these values are identical^ to those derived in reference 10 hy means 
of a different method. 


CONCLUDING REMAE?KS 


The wave drag of any object in a supersonic flow field governed hy 
the wave equation is identical at a fixed Mach number to the average wave 
drag of a series of "equivalent” bodies of revolution. The streamwise 
gradient of area, measured normal to the free stream, at a section of 
each of these equivalent bodies of revolution is given by the sum of two 
quantities: 

1. The streamwise gradient of area, meastured in parallel oblique 

planes tangent to the Mach cones, along the given object 

2. A term proportional to the resultant force on the object measured 

in the same oblique planes. 

Obviously, a formal general application of the above result requires 
a con^ilete knowledge of the shape of the disturbing object and its pres- 
siire distribution - a knowledge which always, of course, fixes the drag 
of the object. However, in many special applications one or the other 
of the two terms mentioned above is smaT T or can be estimated accurately 
enough without a detailed knowledge of the entire airplane or its surface 
pressures. For example, if one wishes to find the wave drag of a wing- 
body combination that is symmetrical about a horizontal, plane (e.g., a 
thin nonlifting wing mounted centrally on a body of revolution) , it is 
not necessary to know the pressures anywhere on the wing, since (if the 
wing is thin enough) their contributions to the resultant forces on 
oblique sections are negligible. Hence, for such an example the force 
term in the drag equation depends only on the pressure over the body and, 

■^In the equation for radial force in Graham's report the exponent 2 
was emitted from the p term. 
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in fact, only on the aBymmetry of these pressures in the ohlique (Mach) 
planes. If the body is slender, the latter effect is small relative to 
the oblique-area gradient, and, for such configurations, the form of the 
supersonic area rule posed by Jones is seen to be a good approximation. 

Ames Aeronautical Laboratory 

National Advisory Ccaumlttee for Aeronautics 
Moffett Field, Calif., Jan. 18 , 1955 
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